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INTRODUCTION 

This Study describes a six - degree -of -freedom 
hybrid simulation for the optimum design of a 
space vehicle re - entry flight control system , 
utilizing the EAI HYDAC® 2400 Combined 
Hybrid Computing System. The specific purpose 
of the simulation is the evaluation of a tempera- 
ture-rate flight control system (TRFCS) which 
utilizes temperature sensors to provide short- 
period stabilization and long-term guidance 
during atmospheric re-entry. 

During the past tenyears, a great deal of research 
and development work has been conducted on 
various types of re-entry vehicles and numerous 
techniques for guiding and controlling such 
vehicles have been proposed. The new flight 
control system described here is, by its nature, 
simple and reliable, and inherently insures a safe 
re - entry. The ensuing discussion explains the 
general re-entry problem and the unique control 
system being studied, and includes a detailed 
discussion of the simulation equipment required 
and its programming. 

GENERAL RE-ENTRY CONSIDERATIONS 

To obtain a better insight into the need for a 
simple yet reliable guidance and control system 
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for lifting re-entry vehicles, abrief discussion of 
the general re-entry problem follows. 

The re-entry vehicle utilized in the study (and 
shown in Figure 1) is an unpowered lifting vehicle 
with wings which, unlike its ballistic counterpart, 
is highly maneuverable and can be landed on 
conventional runways. 

The lift and drag coefficients (Cl and Crj) oi a 
typical high L/D lifting re- entry vehicle are shown 
in Figure 2 as a function of angle of attack. A 
plot of lift to drag ratio (L/D or Cj_/Cd) is also 
shown since it is of importance in determining the 
range of the re-entry vehicle. During the high 
velocity portion of re-entry flight the vehicle 
operates on the high Cl (large a?) side of (L/D) 
max. (15° < oc < 65°). These larger lift coefficients 
yield a re-entry trajectory with lower dynamic 
pressure, acceleration, and temperature. Figure 
3 shows atypical uncontrolled re-entry trajectory 
with a angle of attack, a , of 35 deg. Figure 4 
defines the coordinate system and symbol. 




0( - ANQLC OF ATTACK 

Figure 2. Typical Values of C^, Cp, and L/p for a Lifting 
Re-entry Vehicle 

© Electronics Associates, Inc. 1965 
All Rights Reserved 
1 Bulletin No. ALHC 64041 



During most of the vehicle flight, r G + h « r 0 and 
y * o. Thus, the following approximate equations 
can be used: 
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Figure 3. Typical Re-Entry Trajectory 
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equilibrium glide lines and the boundaries of the 
re-entry corridor. The lower boundary of the 
corridor is determined by the temperature and 
load limits of the vehicle and the upper boundary 
by the recovery ceiling. The recovery ceiling is 
defined as that altitude (with altitude rate, h = 0) 
from which the vehicle can recover without 
exceeding the temperature and load limits of the 
vehicle. 
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Figure 5. Relative Velocity-Thousands of Feet Per Second 
Typical Safe Operating Corridor for Lifting 
Re-Entry Vehicle 



One of the primary re-entry problems is in 
keeping control of the vehicle so that temperature 
and load limits are not exceeded and that a smooth 
equilibrium glide is established. The heavy 
trajectory shown on Figure 5 is a typical uncon- 
trolled re-entry with its familiar skipping oscil- 
lations which cause the vehicle to approach 
dangerously close to the heat limits. The 
temperature rate control system being simulated 
was developed to eliminate these skipping 
oscillations and to reduce the peak temperatures 
during re-entry. 

Another re-entry problem is in carefully managing 
the energy of the vehicle as it re-enters so that 
the desired terminal point is reached. The range 
capability of the re-entry vehicle can be deter- 
mined quite readily in the re-entry portion of 
flight since the vehicle is near equilibrium glide 
where h w 0. If both sides of Equation 2 are 
divided by L/M or g - V 2 /r, the following is 
obtained: 



For a lifting re-entry vehicle, the lift is large 
enough to allow the vehicle to fly along an equil- 
ibrium glide path where h = o. Along this path/, 
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where q is the dynamic pressure (1/2 pV ) and S 
is the wing area. Thus, for a given Cl, there is 
a unique altitude- velocity profile. Figure 5 shows 
the altitude - velocity profile for two different 
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If this is integrated from an initial velocity to 
zero, the range is obtained: 



r o8 



(6) 



Thus, during the re-entry portion of flight, the 
range is a function of velocity and L/D. Range 
control is accomplished, then, by varying L/D, 
i.e.,by varying the angle of attack. From Equation 
6 it should also be noted that the range is very 
sensitive to initial velocity when the velocity is 
nearly equal to orbital velocity, v^g". For the 
re-entry shown in Figure 3, for example, the 
sensitivity of range to initial velocity error is 
approximately 80 NM/fps. 

Lateral maneuverability is obtained by banking 
the re-entry vehicle so that the aerodynamic lift 
vector is rotated, providing a lateral acceleration. 
Figure 6 shows an energy management footprint 
for a typical re-entry flight. The lines of constant 
a. and ft show what attitude must be maintained 
to reach a particular landing site. The dashed 
lines show temperature limits. This large ma- 
neuverability of a lifting re-entry vehicle requires 
a reliable guidance system which will perform 
accurately over the long re-entry, and will 
minimize errors at the desired terminal point. 



The temperature rate control system is being 
simulated to demonstrate its compatibility with 
different types of navigation and guidance systems. 
As will be pointed out in the next sections, the 
TRFCS acts as a filter to the guidance signals to 
insure the safety of the vehicle at all times. 



STATEMENTS OF THE PROBLEM 

Problem Objective: The temperature rate flight 
control system (TRFCS), developed by the AC 
Spark Plug Division of the General Motors Cor- 
poration, is based upon the use of temperature 
sensors instead of conventional inertial instru- 
ments to provide both short-period stabilization 
and long-term guidance during the re- entry flight. 
Details of this flight control system are given in 
the AC Spark Plug Report (1). The mathematical 
formulation for the simulation of the re-entry 
problem was furnished by AC Spark Plug. 

The new control system introduces several 
significant advantages: 

(1) Overall vehicle safety during re-entry . 
This system represents an unorthodox 
approach to the design of an overall 
system of re-entry vehicles. In the 
standard approach, the tendency has 
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Figure 6. Flat Map Showing Typical Energy Management Footprint 
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been towards a complex integrated 
system. In the TRFCS, a successful 
effort has been made to separate safety 
of the vehicle from the task of accurate 
navigation. Because of the inherent 
nature of temperature rate feedback and 
certain selected limits on the control 
authority, the control system minimizes 
skin temperature peaks. The maximum 
"g's" and dynamic pressure are independ- 
ent of initial conditions and maneuvers 
performed. This safety aspect of the 
TRFCS performance is entirely inde- 
pendent of the guidance commands and, 
in fact, the TRFCS serves essentially 
as a filter for them. 

(2) Simple, reliable hardware . This 
separation of control and guidance also 
results in more reliable hardware since 
the failure of the necessarily complex 
guidance system cannot cause the 
complete destruction of the vehicle. 
Furthermore, simple thermocouple 
temperature sensors replace the conven- 
tional gyros and accelerometers. These 
sensors are used to control the flight 
path as well as the short-period 
oscillations in pitch and yaw. The only 
additional sensor required, besides the 
temperature sensors, is a vertical 
reference gyro which, for the safety 
aspects of the re-entry, can be quite 
inaccurate. 

(3) Both manual and automatic modes . In 
case of automatic guidance system 
failure, the TRFCS can be controlled 
manually. The manual flight program to 
be followed by the pilot is very simple, 
and the resulting temperature peaks, 
dynamic pressure, and "g" loads compare 
favorably with those obtained in the fully 
automatic mode. 

During past years, extensive simulation studies 
were conducted by AC Spark Plug to obtain 
familiarity with the control system. A rather 
conventional simulation program was pursued: 
first, analog simulations were performed to gain 
qualitative knowledge of the system, and to 
determine the practibility of this approach; 
second, digital techniques were used to evaluate 
the accuracy of the guidance through TRFCS. 

In the analog simulations, the system character- 
istics were split and analyzed in two independent 
studies. In the first, a three- degree- of -freedom 
simulation of the mass center of the vehicle was 
combined with equations describing the short- 



period pitch dynamics of the vehicle. Pitch axis 
controls and trajectory controls in three dimen- 
sions included an approximate, simple lag 
representation of the lateral response of the 
vehicle. In the second type of simulation, the 
effects of lateral dynamics were obtained by 
simulating the dynamics of the vehicle in detail 
by a standard set of lateral stability equations 
with variable coefficients. The coefficients were 
varied with dynamic pressure, velocity, and 
stagnation temperatures of the vehicle skin, all 
obtained from function generators. The data for 
setting up the function generators came from the 
first type of simulation. In turn, the results from 
the second type of simulation were used to 
determine the lumped lateral response for the 
first type of simulation. Thus, a basis for an 
iterative procedure was established. 

The reason for separating the simulation of the 
pitch dynamics and trajectory control from 
simulation of the lateral dynamics was the 
unavailability of necessary simulation equipment. 
The investigation of aspects of the system such as 
coupling between pitch and roll could not be 
made with the above "split" simulation approach 
and required a complete six-degree-of- freedom 
simulation. 

The next logical step, therefore, was to study the 
system's characteristics in a combined six- 
degree-of-freedom simulation. Here the question 
arose of what computer or computers should be 
used. Past experience had shown that the 
conventional first-analog- then -digital approach 
was definitely not the best. Some of the conclu- 
sions gathered during past simulations are as 
follows: 

(1) Repeatability of analog simulation was 
only marginal (50 miles in range). 

(2) Slowness of digital simulation. Even for 
narrow ranges, determined by previous 
analog simulation, digital simulation was 
too time consuming and, therefore, too 
expensive to optimize parameters. 
Reduction of the digital data also proved 
to be a problem. 

Problem Background: The main objective of the 
problem is to evaluate a TRFCS-controlled re- 
entry both in automatic and manual modes of 
operation. To fully explore the ability of the 
TRFCS to control the short period attitude of the 
vehicle throughout the re-entry, a complete six- 
degree- of -freedom simulation is required. To 
evaluate the ability of the TRFCS to guide the 
vehicle during re-entry to the desired terminal 
point, an accurate and repeatable simulation is 
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required. Economy of analysis should be consid- 
ered, especially in the automatic guidance studies, 
where faster-than-realr-time simulation can be 
employed. 

Computational Requirements: In order to attain 

the above problem objective, the following set of 
rigid computational requirements must be met: 

(1) High accuracy in trajectory calculations 
for the evaluation of the guidance 
capability of the TRFCS. 

(2) Very fast computing capability to 
simulate the high frequency parameters 
faithfully for the short-period dynamics 
of the vehicle. 

(3) Real-time and faster - than- real -time 
simulation for control system evaluation. 
(For economical evaluation of the control 
system in automatic mode, the time scale 
should be as high as possible.) 

On the basis of experience gained during past 
simulations, it was concluded that the simulta- 
neous need for high accuracy and very fast 
computation can only be satisfied by a hybrid 
digital-analog computer. Such a computer would 
allow the programmer to choose either analog or 
digital solution for different portions of the 
problem, trading it with fast processing for high 
resolution, etc. 



PROBLEM MECHANIZATION 



All the advantages of hybrid computation are in 
vain unless very careful consideration is given to 
the programming of the physical system under 
study. This phase of the simulation transforms 
the essentially general purpose computer into a 
simulator of the specific physical system. 



Allocation of Tasks: The first step towards a 
successful hybrid program is the allocation of 
tasks on the computer. The underlying philosophy 
is to subdivide the physical system into sections 
and to assign these to various parts of the 
computer where their speed and accuracy needs 
are best satisfied. As shown in Figure 8, the 
physical system consists of four sections, three 
of them (the vehicle dynamics, the TRFCS, and 
the temperature sensor simulation) constituting 
the attitude control loop, while the vehicle 
dynamics (together with the guidance system and 
TRFCS) form the long period loop. 
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Figure 7. HYDAC 2400 System Configuration 

The HYDAC 2400 Combined System, comprising 
a PACE® 231R-V Analog Computer, a series 350 
Digital Operations System (DOS), and a Series 375 
Digital Computing System, was used to perform 
this simulation. Figure 7 shows a block diagram 
of this system. 

Selection of the HYDAC System enabled the 
programmers to use analog and digital computing 
equipment judiciously and to produce the best 
mechanization of the problem. 
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Figure 8. The Physical System 



The assignment of these sections to various 
elements of the HYDAC 2400 computer is shown 
in Figure 9. (Compare with Figure 7 for task vs. 
function assignments.) The attitude control loop, 
consisting of the vehicle rotational dynamics, the 
TRFCS, and the short period sensor equations, 
are programmed on the analog section. In addition, 
the displays and cockpit simulator are tied into 
the analog since continuous analog signals are 
required. The translational equations of motion, 
long period heat sensor equations, and guidance 
equations are programmed on the 375 because of 
the stringent accuracy requirement. The DOS 350 
provides the master timing, data conversion, 
function generation, and reaction jet control logic. 

The DOS 350 timing and control is essential 
because of the operational differences between the 
analog and digital sections. The analog is a 
parallel continuous computer with computing time 
independent of problem size. The digital is a 
serial, discrete interval computer with computing 
time directly dependent on the size of the problem. 
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Figure 9. Block Diagram of Complete System 



The DOS, through its timing and controls, syn- 
chronizes the calculations on each computer and 
controls the flow of information between sections. 
Function generation and the reaction control jet 
logic are ideally suited to the DOS since these 
operations can be performed rapidly in parallel 
with the general purpose digital computer so that 
the digital computation time is minimized. 



dos 150 Program: Figure 10 shows a detailed 
block diagram of the DOS. The DOS performs the 
five following functions: 1.) Timing, 2.) Mode 
Control, 3.) Data Transfer, 4.) Function Gen- 
eration, and 5.) Reaction jet control logic. These 
functions are described in detail in the following 
sections. 

1. Timing 

Timing is required on a hybrid computer for the 
following reasons: 

(1) To make the mathematical time step used 
in the integration in the digital program 
correspond to the physical time scale 
used on the analog computer; i.e., to 
make the digital and analog sections run 
in synchronism. This synchronism is 
accomplished by sending a periodic 



master time pulse, Tl, which initiates 
the calculations for each time step in the 
digital computer. 
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Figure 10. DOS 350 Block Diagram 
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(2) To time information transfers between 
the analog section and the digital section. 
Not all transfers are at the same rate 
since the serial memories of the DOS are 
used for function generation of aero- 
dynamic moment and force coefficients. 

These variables must be transferred at 
a high rate since the functions are used 
in the short period rotational dynamics 
of the system. On the other hand, those 
variables relating to long term trajectory 
variations are transferred to the digital 
section at alower sampling rate. Timing 
is clearly needed to control these two 
different transfer rates. 
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Figure 11. Master Timer - Block Diagram 



Timing is controlled from the DOS with a Master 
Timer which controls the sequence of events 
occurring per cycle of operation. The master 
timer is a BCD counter driven by the SM8 signals 
occurring at the rate of 2^ pps (approximately one 
pulse every 2 ms.) Figure 11 shows a block 
diagram of the master timer. 

As shown schematically in Figure 12, events are 
scheduled with respect to the master timing 
signal, Ti, (every 25 SM8's), the transfers AD 
and DA being initiated at the specified times TAD 
and TDA respectively. At the end of one compu- 
tation cycle, the digital program jumps to the 
main executive loop and waits there until Ti is 
received. In this manner, the digital section can 
be run synchronously with the analog. For 
instance, if Ti occurs every 50 ms and the time 
scale desired is 2.0 times real time, the time 
step to be used in the digital program is 0.1 sec. 

2. Mode Control 

A very important function of the DOS 350 is to 
control the modes of operation of the system. 
Communication between the DOS and the 375 
occurs through 16 sense lines which are set from 
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the DOS and sensed on the digital section. In the 
other direction, eight flip-flops on the DOS can 
be set from the 375 with special OCP instructions 
(output control pulses). Of these, four can also 
be reset from the digital console. 

All modes are controlled by pushbuttons from the 
DOS 350 and the following list summarizes the 
state and function of the analog and digital sec- 
tions when the indicated pushbutton is depressed. 



(1) IC (initial condition) 

a. 231 R in IC 

b. 375 in IC 

In IC the 375 goes through all computations with the 
exception of the integration routine, 

(2) TIC (type in intital conditions) 

a. 231 R in IC 

h. 375 is ready to accept new initialization data from 
the typewriter. From this mode the program returns 
automatically to the IC loop. 

(3) TYTI (type titles) 

a. 231 R in IC 

b. The 375 types out title block and line headings for 
the 26 variables chosen for print out. 

(4) TRA (transfer only) 

a. 23 1R in IC 

b. 375 in IC 

This mode is for single stepping through the DA and 
AD transfers and was found very useful for problem 
checkout. The Master Timer is stopped and now in- 
cremented only manually by pushbutton action. The 
data is transferred in both directions, one word at a 
time, each time the pushbutton TEST CONV is de- 
pressed. After completion of the DA transfer, the di- 
gital program jumps to the executive waiting loop 
(EXWL) where it waits for another master timing 
signal from the appropriate pushbutton. 



(7) TS (time scale) 

a. 231 R in operate 

b. 375 in operate 

This pushbutton changes time scales, there being 
two arbitrary time scales available, i.e., real-time 
and twenty-times-real-time. 

(8) HOLD 

a. 231 R in hold 

b. 375 in hold waiting loop 

The 375 stays in a waiting loop in the executive pro- 
gram until another mode is selected. 

(9) DUMP 

a. 231 R in hold during actual dump operation, otherwise 
in IC or operate as previously selected. 

b. 375 goes to output routine at periodic intervals de- 
termined by TP, the printout time interval. The digi- 
tal jumps to the output routine and types out the 
present values of the 26 variables. The DUMP com- 
mand can be given in either IC, Hold or OP and the 
system resumes in which ever mode it was at the 
time of execution. 



3. Data Transfer 

The data transfer for this problem is very 
demanding since two basic sampling rates are 
required, one for the short period aerodynamic 
functions and the other for the long period tra- 
jectory variables. The data transfers will be 
discussed in two parts: 1.) the analog-to-digital 
conversion and 2.) the digital- to- analog conver- 
sion. 

a. Analog-to-Digital Conversion . .. . Two variables, a 

and 5 E , are converted every 5 millisec- 
onds since they are used on the DOS for function 
generation. Once every 50 milliseconds, a , /3 , 
and h are converted and transferred to the 375 
for use in the long period trajectory calculations. 
Figure 12 in the previous section shows these two 
different types of conversions. 



(5) DA Test Pattern 

a. 231 R in IC 

b. The 375 goes directly to the DA transfer and back to 
the EX WL, skipping all calculations. A fixed block of 
data consisting of positive and negative maximum 
values (corresponding to ± WOv on the analog) is 
transferred continuously. This was found very con- 
venient for a quick check of the DA conversions. 
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Figure 13. ADC Data Transfer - Block Diagram 



(6) OP (operate) 

a. 231R in operate 

b. 375 in operate 

This is the normal mode of operation of the system. 



A block diagram of the ADC data transfer is shown 
in Figure 13. Upon command from the master 
timing program on the DOS, the ADC control logic 
increments the multiplexer to the proper channel 
and sends a convert signal to the AD converter. 
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When the conversion is complete, the converted 
data is either loaded into a serial memory on the 
DOS for use in the function generation program 
or sent onto the 375 for use in the digital calcula- 
tions. The data transfer from the DOS to the 375 
is accomplished by loading the data into the buffer 
register and setting the parallel input channel 
ready flip - flop on the digital computer . This 
flip-flop is enabled from the DOS by a pulse on 
the TEiline prior to each A-D transfer. The 375 
then inputs the data through its parallel input 
channel. It can be seen that the DOS controls all 
the A-D conversion, thus minimizing digital 
computing time on the 375. 

b. Oigital-to- Analog Conversion . . . Updated values of the 

aerodynamic force and moment coefficients 
which are generated on the DOS are 
transferred to the analog every 5 milliseconds. 
The trajectory variables (altitude, temperature, 
guidance errors, etc.) are transferred from the 
375 to the analog once every 50 milliseconds. 
Figure 12 in the previous section shows the timing 
of these transfers. 

A block diagram of theDAC data transfer is shown 
in Figure 14. The digital -to -analog transfers are 
initiated by the DOS control logic upon command 
from the master timing. If data is to be trans- 
ferred from the 375, the DOS sets the output 
channel ready flip-flop on the 375. This flip-flop 
is enabled from the DOS by a pulse on the TE 0 
line prior to eachD-A transfer. The 375 will then 
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Figure 14. DAC Data Transfer - Block Diagram 

output information to the buffer register and the 
DAC logic on the DOS loads it into the proper DA 
converter. To minimize digital computer time on 
the 375, four data words are loaded into the four 
DAC's (under the control of the DOS). While this 
data is being converted and then demultiplexed on 
the analog under control of the DOS, the 375 is 
formatting the next four words to be transferred. 
This process is repeated until all sixteen words 
are transferred. Using this technique, the total 
processing time consumed for the transfer 
operation (exceptfor formattingjis only 200 ^secs. 
The force andmoment coefficients are transferred 
by loading the four DAC's from the serial 
memories which store the latest computed values 



of these coefficients. The data is then converted 
and demultiplexed on the analog and then the last 
four coefficients loaded. The DOS controls these 
conversions so that they occur at the end of each 
function generation period and therefore do not 
interfere with the conversions from the 375 to 
the analog. 



4. Function Generation 

The present simulation requires the generation of 
eight aerodynamic force and moment coefficients. 
Among these are four functions of one variable 
(C^ (oc), Ci/s (oc), Cvsa (cc), Ci & (*)), and four 
functions of two variables, (Cm(oc, 8E ), Cl(cc, M), 
Ce^", M), Cy /3( oc , M)). oc, M, and §e are angle 
of attack, mach number, and elevator deflection 
respectively. These functions are generated on 
the DOS for the following reasons: 



(1) The functions of two variables are 
extremely difficult to generate in the 
analog section and would, at best, require 
a number of sums and products of 
functions of one variable. The functions 
of one variable could be generated on the 
analog but were programmed on the DOS 
because of the ease of setup and the 
speed at which the functions could be 
changed to study other vehicle config- 
urations. 

(2) The functions are also very difficult to 
generate on the 375 because of the fast 
sampling rate required. The sampling 
rate on the functions should be at least 
10 samples per second since they are 
used in the short period attitude loop of 
the vehicle. In the real time simulation, 
the functions should therefore be 
sampled at least every 100 milliseconds 
and in the twenty-times-real-time mode 
the functions should be sampled at least 
every 5 milliseconds. To meet this 
requirement and still use the 375, the 
375 would have to be interrupted to 
compute these functions a great number 
of times during the major computation 
cycle. (In this problem, the major 
computation cycle is 50 milliseconds.) 
It is estimated that the computation 
time required to compute these functions 
on the 375 is about 4 milliseconds; 
therefore, if the program is stopped 
every 5 milliseconds for function 
generation, only one millisecond of the 
five millisecond interval can be spent on 
the solution of the long period problem. 
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Figure 15. Block Diagram of Two Variable Function Generation 



Hence, the major computation cycle 
would be 5 times longer and twenty- 
times-real-time runs would be impos- 
sible. 

By use of the DOS, the function generation at 
sample rates of 10 per sec per function can be 
accomplished in parallel with the digital 
calculations in the 375. 



configurations, the sign conventions, and the 
combination of jets required to execute various 
commands. To make the reaction jet control 
easier, a pulse width modulation scheme is used 
which makes the moment proportional to the error 
signal. To prevent continuous pulsing of the jets 
even for small errors, a deadzone is built into 
the controller. 



Figure 15 shows a general block diagram of the 
function generation technique. The program for 
this problem utilizes 2-SM8's in series and allows 
32 curves of 16 points each to be stored. Linear 
interpolation between points is utilized and for the 
functions of two variables several curves are used 
with linear interpolation between them. For the 
present problem, the four functions of one variable 
utilize 1 curve each and the remaining four 
functions of two variables are generated with 7 
curves each. 

Each variable is sampled once per 2 SM8 cycle, 
i.e., once every 4 msec. In a 20-times-real- 
time run this would correspond to a sampling 
rate higher than 10 per second. 



5. Reaction Jet Control System 

Since portions of the flight are outside the 
atmosphere in regions where the dynamic 
pressure is too small to make use of aerodynamic 
surfaces for control, a reaction jet system is 
required. Figures 16 and 17 show the block 
diagram of the digital controller, the jet 
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Figure 16. Block Diagram of Digital Controller 




COMMAND CONFIGURATION 



COMMAND 


JETS 


COMMAND 


JETS 


C» P|T CH OP 


3 




6 


C e PITCH DOWN 


4 


Cy, C<J 


2 


YAW RIGHT 


2, 6 


Cy, , C^ 


1 


Cy YAW LEFT 


•,5 


<>. c ; 


5 


CJ ROLL RIGHT 


2,5 






ROLL LEFT 


1,6 







Figure 17. Relationship Between Commands and Jets 

Figures 18 and 19 show the circuitry required in 
the roll axis control to test the magnitude and sign 
of the roll error and to provide the pulse width 
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modulation. Figure 20 shows the logic required 
to decode the commands so that the proper jets 
are used. 
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Figure 18. Circuit for Testing Roll Error Magnitude and Sign 




Figure 19. Pulse Width Command Modulation Circuit for 
Roll Angle 
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Figure 20. Commands-to-Jets Logic 

The logic operations required in the simulation 
of the control system are performed on the DOS 
in parallel to all other operations in the digital 
section. The simulation of such a system by 
conventional analog techniques is a formidable 
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task (dozens of switches and relays would be 
required). The use of the 375 for such an 
operation would result in a significant increase 
in digital computation time because of the large 
number of logic operations and the fast sampling 
rate required. 



Digital Calculation on the 375: The digital program 

was written to make the digital calculation time 
as small as possible and to stay within a 
memory capacity of 4,000 words. The material 
which follows gives a description of the equations 
which are solved on the 375 and the details of the 
digital computer programming. 

1. Summary of the Digital Calculati ons 

Figure 21 shows the block diagram of the system 
of equations to be solved on the 375. A complete 
summary of the equations and symbols is given in 
Appendix A. 




ab'c' -.space fixed reference frame 
abc ;earth fixed reference frame 
ned -.local horizontal reference frame 
x y 2 : body reference frame 
\or : spherical coordinates relative to 
earth fixed reference frame 

Figure 22. Orientation of Space Fixed, Earth Fixed, and Body 
Reference Frames Used in Six Degrees of Freedom 
Simulation 

The three — degree - of - freedom translational 
equations of motion are solved in a local 
horizontal coordinate system with axes along the 
north, east, and radial directions as shown in 
Figure 22. The gravity and altitude calculations 
are based upon an oblate model of the earth, and 
the U.S. Standard Atmosphere (1962) was stored 
in table form on the 375. Most of the equations 
are conventional and quite straightforward, with 
the exception of the heat transfer, temperature, 
and guidance equations which will be discussed in 
more detail. 



The heat transfer and temperature equations are 
as follows: 

2.70893 /V \ 2 

*..-"«■>"©' 

4 C = q^COS^cc 
^R= <5 R S COS60C 

T s=[^ (q RS +qcs)] 1/4 
T = [ ^r + V ] ' /4 

T S = _(3 V+ — ) 



The quantity p e is calculated from the stored 
density table by use of a second order curve fit 
technique to determine the derivative. These 
equations were calculated in the digital section 
because of the large numbers of complex functions 
required (i.e., sine, cosine, log, exponential, 




HEADING ERROR: (T^-CT 

Figure 23. Definitions of A, 0 and Range Angles 
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square roots, etc.). However, they are easily 
computed on the 375 by use of library subroutines. 
Fortunately, these terms were also slow-varying, 
trajectory-dependent quantities so that the slower 
serial computation on the digital is acceptable. 
The temperature rate equations involving angular 
rates are calculated on the analog because of their 
rapidly varying characteristics. 

Because of their critical accuracy requirements, 
the guidance calculations are also performed on 
the digital section. These can be separated into 
two parts: (1) the determination of the range- to- 
go (RTG), heading angle to target {a >p), cross 
range -to -go (CRTG), and heading error and, 
(2) the generation of down range and cross range 
errors (DRg and CRj;) which are sent to the 
TRFCS and subsequently controlled to zero. The 
equations for range-to-go, heading angle, cross- 

range-to-go and heading error, which are 
determined from spherical trigonometric re- 
lations, are given in Appendix A, Figure 23 
illustrates the physical significance of these 
quantitites. The down range and cross range 
errors used for guidance are generated as 
follows: 



DR E = RTG - RTG j 



CR E = 

RTG d 
CRTG 



CRTG - CRTG, 



f,(V) 



The quantities f i (V) and f2 (V) are the values of 
RTG and CRTG for a nominal re-entry. They 
are stored as a function of the relative velocity. 

2. General Description of Digital Program 

In addition to solving the translational, tempera- 
ture, and guidance equations, the digital program 
must accept mode control and timing commands 
from the DOS, scale variables as they are 
transferred in and out of the digital section, and 
finally, provide various input-output functions as 
described in the DOS section. 

Figure 24 shows the flow diagram of the digital 
program. The executive program is discussed in 
detail in the section that follows. The dashed lines 
in the diagram link together the calculations which 
are performed when the IC mode of operation is 
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Figure 24. Flow Diagram of Digital Program 
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selected. All operations are performed except 
the Runge-Kutta integration shown in block 9. 
When the operate mode is selected, the integration 
loop is entered and three passes around 
the loop are made (see description of numerical 
integration which follows). At the end of three 
passes, the positions and velocities, temperatures 
and guidance functions are updated and theanalog- 
to - digital and digital - to - analog transfers are 
made. The executive loop is then re-entered and 
the cycle continued until a different mode is 
selected. The following section discusses this 
executive loop in detail. 



3. Digital Executive Program 

The main medium of communication through which 
the 375 receives the mode commands and timing 
from the DOS 350 is the digital executive program 
shown in Figure 25. This program is essentially 
a chain of test instructions through which the 
digital section interprets the DOS 350 commands 
and executes them by jumping to the respective 
portion of the stored program. 




Figure 25. Digital Executive Program 



Since some of the symbols are not identified in 
the figure, a short explanation of the executive 
program follows. The first test determines 
whether T is greater than TF, where T is the 
current time (time elapsed since the beginning 
of this simulation run) and TF is final time 
indicating the end of the run. If final time is 
reached, results are outputed and the digital 
section is halted. When restarted, the digital 
program begins with the next decision. The next 
decision determines whether fast or slow time 
scaling is requested. If fast time scale is 
selected, the increment size for integration, DLT 
(delta t), becomes 1 sec; otherwise, it is 0.1 sec. 
The subsequent two test instructions simply 
switch the digital program into TIC or TYTI 
modes if so commanded by the mode control 
logic. (See Section B for a definition of these 
modes). The T > TT decision determines if it 
is time for the periodic dump of pre - selected 
digital parameters. The variable TT holds the 
time for the next dump, say 20, 30, 40, etc. 
seconds. TP is the printout time interval and 
increments TT when the dump time is reached. 
If the dump switch is on and T > TT the output 
dump is performed; if not, the program proceeds 
to the next decision. The next test keeps the 
entire digital program under the timing control 
of the master timer. The 375 cannot proceed 
further until the next Tl pulse indicates the 
beginning of the next computational cycle. 
Further testing takes place only after the Tl 
pulse arrives. The TDA test enables the 
programmer to place the digital program into a 
loop where only DA transfer is performed. This 
mode serves as a convenient test for DA 
conversion. (This instruction may be removed 
from the program after the check-out phase is 
over). The last few tests are self-evident and 
require no explanation with the exception of the 
seemingly superfluous second test for dump. The 
test for dump in hold mode makes it possible to 
dump any time by putting the computer in hold 
mode prior to the dump request, in addition to or 
in lieu of periodic dumps. In this manner the 
programmer can determine parameter values 
with digital accuracy at any time during the 
simulation. 

4. Numerical Integration 

Unquestionably, one of the most complex and 
time-consuming parts of the digital program is 
the solution of the six simultaneous differential 
equations which provide the position (altitude, 
latitude, and longitude) and the velocity (radius 
rate, velocity east, velocity north) of the vehicle. 

The numerical technique selected for this 
calculation was the fourth order Runge-Kutta 
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method. The basic method as applied to a single 
differential equation is described briefly. 

Let ^ = f(x,y) represent any first order equation 

and 

h K 2 
K 3 = hf(x n + T , y n + T ) 

K 4 = hf(x„ +h, y n + Kj ) 

A y =-g-(K, + 2K 2 + 2K 3 + K 4 ) 

Then, x n+1 = x n+h and y n+1 = y n+ Ay 



The preceding technique was extended to a system 
of six equations in an obvious manner. 



While preparing the digital program an effort was 
made to combine the calculations due to the 
integration with other necessary calculations to 
minimize the processing time as well as the 
memory space requirements. This approach 
makes it difficult to trace the integration on the 
flow diagram of the digital program shown in 
Figure 24. For example, the Runge-Kutta 
constants are determined in block 6 and the 
variables are incremented for the calculation of 

the next constant in block 9 (such as V E + ^ -> 

VE^etc, and the equations are evaluated with the 
incremented variables in blocks 3, 4, 6, and 7. 
The final calculation of the variables is performed 
in block 9. 



The increment for the second interval is computed 
in a similar manner by means of the same 
formulae. 

A computer flow diagram for this scheme is shown 
in Figure 26. The loop on the left side of this 
diagram corresponds to the constant calculations 
(ki,k2,k3,k4), while the loop on the right updates 
the independent variable and starts the calculation 
of the next increment. When the number of 
increments computed (n) is equal to the number 
desired (nf), the computation ends. Note that x s 
and y s , the starting values for the current 
increment calculation, must not be destroyed 
during this process until the new starting values 
are produced since all formulae are dependent 
upon them. 




The question may be asked by the competent 
reader: Why was this particular scheme chosen 
from many others available? In order to answer 
this question, the selection criteria are listed 
below. 

(1) Self-starting method. It was desirable 
to choose a self - starting method to 
simplify programming and reduce mem- 
ory requirement. (If the method is not 
self- starting, another method is needed 
to calculate the first few points of the 
solution. This virtually doubles the 
integration program.) 

(2) Accuracy. The error of the fourth order 
Runge-Kutta method is of the order of h 5 
and provides a wide enough range for h 
within which the accuracy of the calcula- 
tion is acceptable. (During the simulation 
runs this reasoning proved to be a valid 
one.) The accuracy of calculation using 
the second or third order method was 
considered to be insufficient or marginal. 













EVALUATE 
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GET ki=fh 




Figure 26. Block Diagram of Numerical Integration 
(Fourth Order Runge-Kutta Method) 



(3) Large time steps in integration. While 
processing time associated with this 
method is considerable, the accuracy is 
good enough to allow larger time steps. 
The increase in computation time using 
Runge-Kutta is more than off set by the 
increase in allowable step size. 

(4) Ease of changing step size. In the 
selected method, the length of the step 
canbe modified atany time in the course 
of the computation without additional 
labor. This was considered to be a 
substantial advantage since the simu- 



15 



lation is frequently called upon to change 
time scale, thus requiring a different 
time step for the integration. 



of the equations simulated in the analog section 
is given in Appendix A. The symbols used in the 
equations are defined in Appendix B. 



5. Utility Programs 

In order to satisfy the various data handling needs 
(type in, type out, punch tape, convert binary to 
decimal, etc.), an excessive amount of digital 
programming must be done. Fortunately, all these 
programs are already available, tested, and 
clearly described in the 375' s "software package". 
This is no small feat if one considers that the 
above programs, together with the numerous 
subroutines, normally add up to about 75% of all 
digital programming required. Furthermore, as 
it will be pointed out later, several programs are 
available to make the debugging and updating 
procedures efficient and fast. 

Analog section: The analog computer is the one 
link of the HYDAC 2400 system ideally suited for 
control system simulation by virtue of its cap- 
ability for high speed and parallel computation, 
and its input/ output flexibility. Output data can be 
displayed in a multitude of forms such as X-Y 
plots, strip chart plots, oscilloscope displays, 
auxiliary meters, etc. Special purpose input 
equipment can be adapted easily for compatible 
operation with the analog computer. 

1. General Mechanization 

The block diagram of Figure 27 delineates the 
mechanization of the analog program. A listing 
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As previously mentioned, the TRFCS and rotational 
dynamics were programmed on the analog due to 
their rapidly varying characteristics. When the 
vehicle is in the atmosphere, the vehicle attitude 
is controlled by aerodynamic control surfaces. 
The aerodynamic control moments are calculated 
from the surface deflections and the moment 
coefficients generated on the DOS 350. Out of the 
atmosphere, the vehicle attitude is controlled by 
areaction jet control system. The jet pulse logic 
is generated on the DOS 350 and transferred to 
the analog section where the moments are 
produced and fed into the angular acceleration 
equations. The angular rates, generated from the 
angular acceleration equations, are used to 
calculate the Euler angles 0 and J0. These angles 
are used to resolve gravity into the body axis for 
use in the a and /3 equations shown below. 



g V, L 
oc = {-—— — ) (cos 6 cos d> cosoc+ sin(9 since)— — 
V rV mv 



+ Q — (Pcos qc+ R sin cc)/3 



gr V, 

_D 

f . mV 



!= — — — (cos 6 sin c6) y- Rcosoc+ Psinoc+ — y 

vi 2 

— (gr — (cos 9 cos <f> sin oc— sin 6 cos oc) /3 



These seemingly redundant force equations, cc and 
/3, are computed on the analog since they are 
required in the short period sensor equations 
shown below. The aerodynamic force coefficients 
in the i and /3 equations are generated on the 
DOS 350. The geometric definitions of <* and p 
are shown in Figure 28. 




xijz — body Axes 



Figure 28. Definition of « and (3 Angles 



The short period sensor equations are: 



Figure 27. Block Diagram of Analog Section 



T = T s (l - 1875 cc ) - .375ccccT s 
AT = .021(T s /3 + T s jS) 
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The quantity T is the temperature rate at the nose 
sensor; AT is the temperature rate differential 
between the two wing sensors (See Figure 1). The 
stagnation point temperature information, T s and 
T s , which is used in these sensor equations, is 
calculated in the digital section. The T and AT 
are used in the TRFCS control equations shown 
below: 

5 E = f,(T) 

8 a '= K t(fi-/" C ) + K 2 P + K 3 5 R 
Sr = f 2 (AT)+ K 4 5 a 
= f 3 (T) 

The T and AT terms supply damping to the control 
equations, alleviating the heating problems 
associated with an undamped trajectory. Closed 
loop guidance is achieved by adjusting the pitch 
axis controls with a compensated down range 
error, DRg , and adjusting the n c to compensate 
for the cross range error, CRg. The pilot can 
control the range manually by adjusting his 
temperature rate profile to eliminate the displayed 
down range and cross range errors. 

2. Cockpit Simulator and Display Equipment 

A cockpit simulator is utilized to evaluate the 
TRFCS in the manual mode and is trunked directly 
to the analog computer. Computer outputs drive 
display meters on the TRFCS CONTROL PANEL* 
which monitor the following parameters: 

(1) temperatures 

(2) temperature rate 

(3) velocity 

(4) acceleration 

(5) pitch damper 

(6) yaw damper 

(7) bank angle 

(8) heading angle 

(9) elevator trim 

(10) range 

(11) range error 

(12) cross range 

(13) cross range error 

(14) altitude 



*Refer to Figures 28 and 29 of Reference (1). 



As seen in the analog program block diagram, 
many of these parameters are transferred from 
the digital section. Some of the parameters, 
though not necessary for control, indicate tra- 
jectory status and, therefore, maintain the pilot's 
confidence in his control information. The above 
parameters and other pertinent data are recorded 
on strip charts and X-Y plots for permanent 
record of each flight. 



Vehicle attitudes commanded by the flight control 
system are displayed on a large oscilloscope. An 
illustration of the display at various attitudes is 
shown in Figure 29. 




ANGLE OF ATTACK =0 BANK ANGLE =0 




ANGLE OF ATTACK = 45° BANK ANGLE=0 




ANGLE OF ATTACK=45* BANK ANGLE=30° 

Figure 29. Oscilloscope Display of Re-Entry Vehicle 
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APPENDIX A 



I. Digital Computer Equations 

A. Basic Translational Equations in Local 
Horizontal System 

The translational equations are solved in a local 
horizontal coordinate system with axes along the 
north, east, and radial directions. Figure 23 shows 
this coordinate system. 



r = — -g r + (£cosO) V E + 0V N (1) 



D. Velocity, Velocity Heading and Flight Path 
Angle Calculations 

The total relative (to the earth) velocity is calculated 
from the velocity components in the N, E and radial 
directions. The horizontal relative velocity has two 
components, east and north. The easterly relative 
velocity is the difference between Vg and the tangential 
rate of the earths rotation. The relative northerly 
velocity is the same as the northerly inertial velocity. 
The flight path angle, /, is the angle between the 
horizontal and the total velocity vectors. The velocity 
heading angle, a , is measured from east to north. 



V N = — + g N - sinQ)V E -Q; (2) 



V E = _ + sin 0) V N - (<f cos S2) r (3) 



B. Aerodynamic Force Calculation and 
Transformation 

The aerodynamic forces are calculated with respect 
to the relative velocity. vector (lift, drag, and y axis 
acceleration) and then transformed into the local 
horizontal coordinate system. 



V = J V 2 + v 2 

H » E R N I 



V ER = V £ - ro> 0 cosfi 



sintfK =V 



> f = V 



jf= tan 



V 

-1 1_ 



(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



L = g SC L 
D = gSC D 
Y = gSC y 



F L Y D . 

— - =( — cos/it- — sin fj.) cos — smf 
m m m m 



(4) 
(5) 
(6) 

(7) 



..il 



(20) 
(21) 

(22) 



f n , L ■ Y ^ u L Y • ^ ■ / 

— - = -(. — sitifi -t — cos^tj cosct- |A — cosfi - — sin/xj sin / 

m m m mm 

H cos r j smff 



F E ,L . Y L Y . 

— = \ — sin/i + — cos/n) smcr - \\— cosfi sin/i) smy 

m m m mm 



H COS**] COSCT 



(8) 



(9) 



C. Angular Rate and Position Calculation 

The derivatives of the space fixed coordinates, 9. 
and $ , are solved as functions of the inertial 
velocities, V N and V E . The derivative of the 
longitude is calculated by subtracting the earth's 
rotational rate from £ . Since the earth rotates in 
the longitudinal direction, is the latitude. Altitude 
is calculated with respect to an oblate earth. 



12= — 



h = r - r E + 2(1 -e 2 ) sin2{1 



(10) 

(11) 
(12) 

(13) 



E. Density and Dynamic Pressure 

The density is tabulated in square root form to com- 
press the range of variation and was taken from U.S. 
Standard Atmosphere (1962). The dynamic pressure 
is a function of density and velocity. 



vT"= v7( h , 



g = v) 2 

V 2 



(23) 
(24) 



F. Gravity Calculation 



The gravity components are calculated for an oblate- 
spheroid-shaped earth. Since the earth is assumed to 
be rotationally symmetrical, g E = 0. 

.2 

r 2 3 



GM 

gR= 5" L 1 



3Jl E (I - sin 2 ft)] 



2GMJr g . 

-. sinU cosU 



(25) 
(26) 



G. Heat Transfer and Temperature Equations 



These equations are generated in the digital program- 
and transferred to the analog for calculating para- 
meters in the flight control equations. 
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2.70893 
8cs = — — VT 



g RS = 144. 9Rp 



8c = gcs^ os 



g = g cos 
R RS 



T s = t^- (g RS + .g cs ) ] 



T = [— (g R + g c )l 



„ T S V 0 i 

g T S= + g cs ) dt 

Pe p dn 



(27) 

(28) 

(29) 
(30) 

(31) 

(32) 

(33) 
(34) 
(35) 



H. Target Heading Angle, Range-to-go and 
Cross Range-to-go 

The range and heading angles are calculated using 
spherical trigonometry and the resulting angles are 
multiplied by the earth radius yielding the range. 

cosO tanflj - sin 0 cos(\ T - a) 
I sin(A T - A) 

sin (A -A) 



(36) 



tan 6 



sin (A T - A) 



RTG cos<7 T [tanfl T sinfl + cosfl cos(Aj-A) 



RTG = 



^RTG^ r E 



tanCR=- tan(cr T -<j) sin0 RTG 
CRTG = (CR) r_ 



I. Heading Error 

The angle between the great circle heading to the 
target, o" T , and the velocity heading, <r . 

Heading Error = a - <j 

J. Guidance Equations 

Error equations for use in the TRFCS 



DR E = RTG - RTG d 

RTG d = f, (V) 

CR E = CRTG - CRTG d 
CRTG d = f 2 (V) 



(37) 



(38) 



(39) 



(40) 



(41) 



II. Analog Computer Equations 

A. Angular Acceleration 

The following equations are calculated in the body 
system: 



M ■ 



X 

M+ M, 



fly - I X \ N + 



(42) 
(43) 

(44) 



B. Aerodynamic Torques 

These parameters are generated by the moment of 
aerodynamic forces. 

L=gSb[C |jQ ( K )^ + C l5r 8 r+Wa ] 



(45) 
(46) 



M= g Sc [C m («,fi E )] 
N = g Sb [C^Cac) B + Cr, Sa (oc) S + C^ (oc) 8J (47) 

C. Euler Angles 

Generated for local horizontal to body axis 
transformation. 



<9 = /(Qcosc/>- R sia(p) dt 

0= /[P + tan0(Qsin<£ + Rcos</>)]dt 

D. Attitude Rate Equations 

Generated in this form for use in the heat sensor 
equations. 



« = .X!_Ycos 6 cos (bcosa + sin 6 sinoc) - ^ 



(48) 
(49) 



mV 



+ Q - (Pcos a + R sin oc)B 



(50) 



(52) 



cos d sin eh ) - R cos oc +P sin cc+ — — 

mv 

2 (51) 

+ — B - (e. (cos 6 cos eh since- sin 6 cos oc) B 

mv r 

Ji = P cos « + QB + R sin « 

E. Heat Sensor Equations 

These short period equations are used for control 
and are discussed in the text. 

f = T s (1 - .1875oc 2 ) - .375 oc ce T s (53) 

AT=.021T sl 8 (54) 

AT= .021(T S /3 + T s /3) (55) 
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F. TRFCS Control Equations 

The control parameters include the variable surfaces 
and the commanded bank angle. 



S E = f^T) 



(56) 



8 a =K 1 ( / z-, lc ) + K 2 P+ K 3 S R 
5 R =f 2 (AT) + K 4 S a 



(57) 
(58) 
(59) 
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APPENDIX B 
DEFINITION OF SYMBOLS 



Variable 



Definition 



Variable 



Definition 



M 

a 

i 

A 

<x 

k 
P 

Q 



acceleration north (local level 
coordinates) 

acceleration east (local level 
coordinates) 

acceleration up (local level coordinates) 



aerodynamic acceleration north 

aerodynamic acceleration east 

aerodynamic acceleration up 
gravitational acceleration up 
gravitational acceleration north 

lift acceleration 

drag acceleration 

side force acceleration 

velocity north, inertial 

velocity east, inertial 

radius rate 

total inertial velocity 
total relative velocity 
relative horizontal velocity 



velocity north, relative (equals V ) 
velocity east, relative 



mach number 

latitude rate 

longitude rate, inertial 

longitude rate, relative 

rate of change of angle of attack 

rate of change of angle of side slip 

roll angular acceleration 

pitch angular acceleration 



h 

DR E 
CR 

E 

RTG 



yaw angular acceleration 
roll rate about the body x-axis 
pitch rate about the body y-axis 
yaw rate about the body z-axis 

Euler rates of body axis with 
respect to local horizontal 

aerodynamic torque about the body 
x axis 

aerodynamic torque about the body 
y axis 

aerodynamic torque about the body 
z axis 

reaction jet moment about the body 
x axis 

reaction jet moment about the body 
y axis 

reaction jet moment about the body 

z axis 

latitude 

longitude 

angle of attack 

angle of side slip 

bank angle (lift vector with vertical 
plane) 

body roll angle 
body pitch angle 

body heading angle 
bank angular rate 
elevator deflection angle 
aileron deflection angle 

rudder deflection angle 
flight path angle (+ for r +) 

velocity heading angle (to the north 
from east) 

heading to target (+ to the north from 
east) 

air density 

radial distance from earth center to 
vehicle 

altitude 

down range error 

cross range error 
range -to -go 
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Variable 



Definition 



Variable 



Definition 



CRTG 
C 



1/8 

\ 
\ 
\ 



cr os s -r ange -to -go 
side force coefficient 

dynamic pressure 
lift coefficient 

drag coefficient 

pitching moment coefficient 

rolling moment coefficient due to 
side slip 

yawing moment coefficient due to 
side slip 

rolling moment coefficient due to 
aileron deflection 

yawing moment coefficient due to 

aileron deflection 

yawing moment coefficient due to 

rudder deflection 

rolling moment coefficient due to 

rudder deflection 

radiative heat transfer rate at 
stagnation point 

convective heat transfer rate at 
stagnation point 

radiative heat transfer rate at sensor 
location 

convective heat transfer rate at sensor 
location 

total integrated heat flux at stagnation 
point 

temperature rate at sensor location 
temperature rate at stagnation point 



T 
T 

s 

AT 

GM 
J 



temperature at sensor 
temperature at stagnation point 

differential temperature rate 
between wing sensors 

earth gravitational constant 

earth oblateness factor 

radius of the earth-equatorial 

eccentricity of the earth 
rotation rate of the earth 

emissivity 

Stefan Boltzman constant 
sea level gravity 
mass of the vehicle 

wing area 

wing span 
wing chord 

moment of inertia about x axis 
moment of inertia about y axis 

moment of inertial about z axis 

nose cap radius 

exponential index of the atmosphere 
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